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Lattice results relevant to the CKM matrix determination 

Damir Becirevic®' 

^'Dipartimento di Fisica, Universita di Roma “La Sapienza”, Piazzale Aldo Moro 2, 1-00185 Rome, Italy 

A brief discussion of the recent lattice studies of the quantities that are relevant to the Standard Model 
description of CP-violation in the hadronic sector is presented. A comment on 5 —> A *7 decay is made too. 


1. Introduction 

The most popular way to estimate the amount 
of CP-violation in the Standard Model is to study 
the shape of the CKM unitarity triangle involv¬ 
ing the third generation of quark flavors. The 
standard unitarity triangle analysis (UTA) [Q, 
which requires the least theoretical assumption, 
comprises the and B^ — mixing am¬ 

plitudes, as well as the SU(3) breaking ratio in 
the B^ — B^ system. Additional constraint comes 
from information on |14b/t4b| ratio, which will 
be discussed towards the end of this review. The 
three “golden” mixing amplitudes over-constrain 
the CKM-triangle through a hyperbola and cir¬ 
cles in the {p, 77)-plane which emerge from the 
comparison of the theoretical predictions with the 
experimental determinations of |eiy|, Anid and 
the limit on Amg/Amd- Besides the CKM pa¬ 
rameters, the theoretical expressions contain the 
Wilson coefficients and the hadronic matrix el¬ 
ements. The Wilson coefficients encode the in¬ 
formation about physics at high energy scales, 
they are computed in perturbation theory and 
are known at NLO accuracy (for a review see 
e.g. iQ). The hadronic matrix elements, instead, 
describe the low energy QCD dynamics and their 
computation requires a theoretical control of non- 
perturbative QCD. To be specific, the following 
matrix elements need to be computed: 

{K°\Ogd{f^)\K°) = 

( 1 ) 

{B°\OM\B°g) = ^mljlBBAf^), 


where the operator Oqq = {Qq)v-A {Qq)v-A- 
On the r.h.s. of eq. (|^ we recognize the 
pseudoscalar decay constants /k and , 

and the familiar “bag” parameters, Bk{^j) and 
That is where the lattice QCD enters 
the stage and provides us with a completely non- 
perturbative way to calculate these (and other) 
hadronic quantities. The following two points 
should be emphasized: Q 

(g) Lattice QCD is not a model but the ap¬ 
proach allowing the computation of vari¬ 
ous Green functions from hrst principles 
of QCD: the only parameters entering the 
computation are those that appear in the 
QCD lagrangian, namely the quark masses 
and the (bare) strong coupling; 

(g QCD is solved numerically on the lattice 
by using the Monte-Carlo methods which 
induce the statistical errors in the results. 
The central limit theorem applies and en¬ 
sures that those errors fall as l/VA’conf, 
with A^conf being a number of independent 
SU(3) gauge held conhgurations on which 
a given Green function has been calcu¬ 
lated. Thus, in principle, computation of 
the hadronic quantities on the lattice can 
be made to an arbitrary accuracy. 

In practice, however, various approximations 
have to be made, each introducing some system¬ 
atic uncertainty into the hnal results. Yet all 
those approximations are improvable by: (a) in¬ 
creasing the computing power so that, one by one, 

^For details about formulating QCD on the lattice, please 
see e.g. g. 
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all the approximations can be removed from the 
computation; (b) improving the approach the¬ 
oretically so that the hadronic quantities, com¬ 
puted on the lattice, converge faster to their val¬ 
ues in the continuum limit. 

Now I return to the lattice determination of the 
“golden” hadronic quantities entering the stan¬ 
dard UTA. The present situation is discussed at 
length in ref. [Q, and the actual results are: 

.Bk = 0.86(6)(14), 

Bb, = 1.34(12), 

/b,= 203(27) MeV, (2) 

fBjfB, = 1.18(4) (til) , 

e = (/s. V^)/(/b,V^) = 1.18(4) . 

The above quantities are completely dominated 
by the systematic uncertainties, some of which 
are discussed in what follows. 

2. — AfO mixing 

Results of the precision lattice computations 
are presented as B^^{2 GeV), in the quenched 
approximation, i.e. the quark loops in the back¬ 
ground gauge field configurations are being ne¬ 
glected, n/ = 0. In the real physical world, how¬ 
ever, one needs n/ = 3. Numerically, the conver¬ 
sion from the MS(NDR) scheme to the renormal¬ 
ization group invariant (RGI) form, Bk, is almost 
completely independent of n/ (see discussion in 
ref. |]). By using = 338(40) MeV §, I 

obtain 

Bk = 1.382(15) . B^{2 GeV). (3) 

A map of the precision quenched lattice QCD de¬ 
terminations of Bx is shown above the dashed 
line in fig. |^, where the asterisks denote the re¬ 
sults presented this year. The most accurate is 
the value obtained by using the so-called stag¬ 
gered fermions j|], 13x — 0.87(6). The advan¬ 
tage of that action is that it preserves the chi¬ 
ral symmetry, but the drawback is that it breaks 
the flavor symmetry. With Wilson fermions, in¬ 
stead, the flavor symmetry is preserved but the 

^For a comparison, with ^ = 250(25) MeV, one has 
Bk = 1.388(8) . 5^(2 GeV). 
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Figure 1. Quenched lattice results for the Bk parameter 
from 4 different lattice actions. The values above/below 
the dashed horizontal line are obtained after/without ex¬ 
trapolating to the continuum limit. The second result 
by SPQcdR, is obtained by using the method proposed 
in ref. |ll[ . 


chiral symmetry is sacrificed. The loss of the chi¬ 
rality induces the additive renormalization of the 
operator Oad^ which is more involved than the 
(usual) multiplicative one. The non-perturbative 
method to subtract these mixing has been first 
implemented in ref. by using the quark Ward 
identities, and then in ref^0 by using the so- 
called RI/MOM method [|1()[. The problem of 
spurious mixing has been alleviated by a judi¬ 
cious application of the Ward identity on the par¬ 
ity violating operator which does not suffer from 
the spurious mixing problem on the lattice |Tl[] . 
That proposal has been implemented in the high 
statistics lattice simulation and the result, in the 
continuum limit, reads Bx = 0.88(11). The re¬ 
sults below the horizontal line in fig. do not 
refer to the continuum limit. Their major signif¬ 
icance lies in the fact that they are obtained by 
using the actions that satisfy both the chiral and 
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the flavor symmetry. Computations with these 
actions are much more involved and require huge 
computational resources. This year, first results 
for Bk, obtained by using the so-called overlap 
fermions, have been reported by the MILC 
and by the Boston-Marseille jl^ collaborations. 
Even though it is early to draw conclusions, the 
observed good agreement with the values above 
the horizontal line in fig. is encouraging. Some¬ 
what lower are the values obtained by using the 
so-called domain wall fermions (DWF). Although 
their statistical errors are very small, in my opin¬ 
ion, their systematic errors are not fully realistic 
(see discussion in ref. Q). Nevertheless, it should 
be stressed that the results for Bk obtained by 
using DWF agree -within the error bars- with 
other lattice determinations. The remaining sys- 
tematics is almost entirely due to quenching. Al¬ 
though the first (low statistics) lattice studies 
suggest that the unquenching is not changing the 
value of Bk we follow ref. 0 and assume 
that the quenching error is of the order of 15%. 
Thus it is of upmost importance, and challenging 
task for the lattice community, to unquench Bk- 

^s{d) - ^°sid) mixing 

In this case, besides the light quark, one also 
has to deal with the heavy quark on the lattice. 
Since the Compton wavelength of the realistic 
6-quark is smaller than the presently attainable 
lattice resolution, one has to adopt some other 
strategy in order to compute the properties of the 
heavy-light hadrons. Four ways to deal with that 
problem have been implemented in the quenched 
computation of the decay constant ’ and/or 
the bag-parameter BB^^^y 

© Lattice QCD: one can afford to work with 
propagating heavy quarks of masses around the 
physical charm quark mass and then one has to 
extrapolate to the 6-quark by means of standard 
heavy quark scaling laws 0| ; 

© Heavy quark effective theory (HQET) on 
the lattice is very difficult to extend beyond 
the static limit (m& ^ oo). The long stand¬ 
ing problem to renormalize the axial current 
non-perturbatively has been solved recently (see 
eg. @)- The static results will therefore be pos¬ 


sible to use in constraining the extrapolations of 
the results obtained in the full theory, i.e. with 
propagating heavy quarks; 

© Non-relativistic QCD (NRQCD), is the ap¬ 
proach in which the l/rub-corrections are in¬ 
cluded in both the lagrangian and the operators. 
Besides difficulties in renormalizing the operators 
in such a theory on the lattice, a particularly 
problematic is the fact that the expansion is in 
l/(amf,) (a being the lattice spacing), so that the 
continuum limit (a —> 0) cannot be taken; 

© Fermilab (FNAL) approach consists in us¬ 
ing the lattice QCD with propagating quarks and 
“sending” the heavy quark mass over the lattice 
cut-off, followed by a redefinition of the hadron 
masses and reinterpretation of the theory in terms 
of a 1/m^-expansion. In some cases the approach 
is believed to be plagued by “renormalon shad¬ 
ows” 

None of the enumerated approaches is fully sat¬ 
isfactory on its own and, at present, they all 
should be used to check the consistency of the 
obtained results. For example, the results for /b, 
as obtained by using various approaches, agree 
quite impressively (see ref. @). Nevertheless, 
the systematics {within the quenched approxima¬ 
tion) cannot be further reduced until the fully 
QCD based lattice method is devised to compute 
the heavy-light quantities directly on the lattice. 
Otherwise, one should wait for the next genera¬ 
tion of parallel computers to do the job. A candi¬ 
date method for the precision computation of /b 
has been proposed this year in ref. [p0[ , in which 
the step scaling function has been employed to 
make a controllable l/(Lmb)-expansion, with L 
being the side of the lattice box. From an ex¬ 
ploratory study, /b = 170(11)(23) MeV has been 
quoted. The authors promise to reduce the error 
bars very soon and to respond to the questions 
raised in ref. pl| |. 

The effect of unquenching has been extensively 
studied this year by the MILC collaboration . 
They conclude that /b^~^ is about 15% larger 
than /b^~ . Moreover, in their preliminary study 
with nj = 3, they show that the ratio Jb^I fBd 
remains essentially unchanged when one switches 
from n/ = 0^2^3 |^. That brings me 
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to the important estimate of the hadronic pa¬ 
rameter ^ (see eq. (H)). The SU(3) light flavor 
breaking effect in the bag parameter is equal to 
BbJBb^ = 1.00(4) and the whole effect of 

^ being different from unity comes from Jb^ /• 
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Figure 2. Quenched and unquenched (nf = 2) lattice es¬ 
timates for the SU(3) breaking ratio /bs /as obtained 
by three differentjnethods indicated on t hp left. For more 
details see refs. Q (MILC) and (CP-PACS). No 

shift due to the chiral log term is displayed. 


3.1. Chiral logs in /bJ/b^ 

Even if one gets a good handle over the heavy 
quark, there is still an extrapolation in the light 
quark mass to be made. With respect to the 
physical strange quark mass, the light quarks that 
are directly accessible from the lattice simulations 
lie in the range 1/2 < niq/ms ^ 3/2, in which 
one observes that the heavy-light decay constant 


changes nearly linearly under the variation of the 
light quark mass. In other words, is accessed 
directly, whereas an extrapolation in the light 
quark mass is necessary in order to reach Jb^- 
It has been argued recently that if in the stan¬ 
dard (linear) extrapolation one also includes the 
chiral logarithmic terms, the SU (3) breaking ratio 
may get shifted from ^bJJb^ = 1-16, by -1-0.16, 
which is a 100% error on the net SU(3) break¬ 
ing 1^. In ref. the shift can be as large as 
-1-0.24. There is an ambiguity while implementing 
the chiral logs in extrapolation. It is not clear at 
what mass scale the chiral logs dominate over the 
other, higher order, terms in the chiral expansion. 
If one sets it for the light meson masses around 
600 MeV, the extrapolated value will be strongly 
shifted upwards (w.r.t. the result of linear ex¬ 
trapolation) The shift is roughly halved if 

the chiral logs are introduced for the light pseu¬ 
doscalar mesons ^ 350 MeV. To avoid that ambi¬ 
guity, ref. proposes to study the double ratio 

R = 

in which almost all the chiral log corrections can¬ 
cel, and one can make the linear extrapolations 
with a small error. From the available lattice data 
with Uf = 2, the NLO chiral correction in the ra¬ 
tio (|^) seems to cancel too. From that consider¬ 
ation, ref. |2^ concludes that /bJ/b^ ^ fK/U, 
where the latter ratio is known very well from 
the experiments, Jk/!-k = 1.22(1). Therefore, 
the shift is indeed present but it is not as spec¬ 
tacular as previously thought. Similar conclusion 
has been reached by the MILC collaboration . 
They showed that if the chiral logarithms are con¬ 
sistently included in the entire analysis of the lat¬ 
tice data, the shift in fBs/fsd is much smaller, 
namely -1-0.04. 

Notice that the second errors in JbJ fBd 
and quoted in eq. (^) (see also ref. |^), reflect 
a guesstimate of the impact of the chiral logs on 
extrapolation in the light quark mass. As we just 
explained, those errors cannot be too large, and 
I personally believe that through the strategies 
such as those proposed in refs. ^ and |^, the 
lattice community will be able to cut down on 
those errors very soon. 


IBs/ fSd 
fK/U ’ 
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4. \Vcb\ from B —s- D^*'>iv 

The extraction of \Vcb\ from the exclusive 
semileptonic modes, B requires the 

computation of the corresponding form factors. 
The cleanest method for their extraction has 
been proposed and implemented in ref. p^ . It 
has been discussed at this conference by J. Si¬ 
mone , and the reader is referred to those ar¬ 
ticles for more details. It would be very nice if 
the other lattice approaches adopt that strategy 
and compute the relevant form factors through 
the double ratios described in ref. [p^. 


5. |I4b| from B 


In recent years several new lattice studies of the 
B ^ TT form factors appeared pO| , pl|J^ , |^ , p^ . 
The present situation with lattice results is de¬ 
picted in fig. where we plot the two form factors 
that parametrize the hadronic matrix element 


{Tr{p')\Vf,\B{p)) = c+F+{q^) + c°Fo(g2), (5) 


with c_|_/o being the known kinematic factors. 
The range of the transfer momenta available from 
this decay is huge: 0 < < {ms — rriTr)'^ = 

26.4 GeV^. The lattice can be used only for 
large g^’s (small recoil), and therefore the as¬ 
sumptions are necessary if we are to cover the 
entire g^-range. In ref. a simple parametriza- 
tion has been proposed which includes various 
symmetry constraints onto the shapes of the form 
factors [ pTlJSSlI , as well as the kinematic condition 
^+(0) = Tb(0). It is often referred to as BK and 
it reads 


F+{q^) 

Fo{q^) 


c {1 — a) 


(1 - gVm|.)(l - aq'^lm%,) ’ 


(1 - a) 


1 - g2/(/3TO|j.) ’ 


( 6 ) 


thus consisting of three parameters only: c, a and 
(3. Results of the ht of the lattice data to this form 
are given in table |^. We notice a pleasant agree¬ 
ment with the results obtained by LCSR p5| , 
which is believed to be the method of choice in 
the low g^-region. We also see a good consistency 
of the results obtained by using three different 
approaches. The agreement is especially good for 
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Figure 3. Recent (quenched) lattice results for the two 
form factors are plotted head-to-head, since they satisfy 
the condition F+(0) = i^o(O)- For illustration a fit to the 
BK-form is shown as well as the prediction obtained 
by using the light cone QCD sum rules (LCSR) ||35| . 

the form factor F+(g^), which is the only one en¬ 
tering the decay rate (for i = e,p.). A bit less 
good is the agreement for the form factor Fo(g^). 
It is worth stressing that so far all the computa¬ 
tions are performed in the quenched approxima¬ 
tion. Recently, in ref. it has been pointed 
out that the quenching effects may be significant. 


Table 1 

Results by various lattice groups are fitted to the form (^) . 
For comparison, the LCSR results are presented in the 
same form (last line). Note that instead of the parameter 
c, we quote F{0) = c (1 — a). 


F{0) a (3 

ref. 

0.30(ti)(th 0.46(t°9)(+37) 

0.28(6)(5) 0.45(17)(;“) 1.20(13)(tJ5) 

0.26(5)(4) 0.40(15)(9) 1.22(14)(tJ5) 

0.33(1^) 0.34(1^) i.3I(t|) 

0.23(14) 0.58(142) 1.28(142) 

UKQCD @ 
APE-i Hi 
APE-II H 
FNAL H 
JLQCD Pi 

0.28(5) 0.32(124) 

KRWWY P 
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particularly for the form factor FQ{q^). That con¬ 
clusion has been reached by confronting the cor¬ 
responding expressions derived in the quenched 
and in the full ChPT. The remaining step, in my 
opinion, is to unquench the B ^ tt form factors. 

In the case of i? ^ piiy, instead of 2 one has 4 
form factors and less kinematic constraints, which 
makes the computation more complicated. Like 
in the case of i? —> tt transition, the lattice re¬ 
sults for B > p decay are accessible for relatively 
large values of q^. The constrained extrapola¬ 
tions to the low g^’s, guided by the symmetry re¬ 
lations |37|j3§|] , are however much more difficult to 
control in this case. To make the phenomenolog¬ 
ically relevant predictions, the lattice results for 

> 10 GeV^ of ref. Q have been interpolated 
and then combined with the light cone QCD sum 
rule results of ref. Q , which are again expected 
to be reliable in the region of low g^’s. From that 
exercise and by using the experimental branch¬ 
ing ratio recently measured by 

CLEO, BaBar and Belle [|l), we obtain 

IKbl = 0.0034(6). (7) 

The reader interested in new results for B —> piv 
decay form factors is invited to consult ref. |40|| . 

6. B ^ K*j 

6 ^ S 7 is a fascinating decay mode: it is me¬ 
diated by the flavor changing neutral current and 
therefore completely dominated by the loop ef¬ 
fects. The hope is that by confronting theory 
Vs. experiment, one can probe the content in 
the loops and perhaps detect the non-Standard 
Model physics contributions. For that purpose 
the accurate Standard Model predictions are nec¬ 
essary. The NLO theoretical prediction for the 
inclusive B —> Xsj decay has been recently com¬ 
pleted in ref. (see also references therein). For 
the exclusive decays, instead, the situation is ob¬ 
scured by the necessity for the accurate compu¬ 
tation of the hadronic form factors 

{K*{p')\sa^^il + j5)b\B{p)) , 

i=l 

where the kinematic factors are known (see 
e.g. H^). When the photon is on-shell, i.e. at 


BHS (1994) 
ref. [43] 

UKQCD-I (1994) 
ref. [44] 

APE (1995) 
ref. [45] 

LANE (1995) 
ref. [46] 

UKQCD-11 (1998) 
ref [47] 

SPQcdR (2002) 

ref [48] 

P.Colangelo et al. 
ref [49] 

P.Ball, V.M.Braun, 
ref [50] 


0 0.1 0.2 0.3 0.4 0.5 



T(0) [B^Ky] 


Figure 4. Lattice results for the form factors (0). 

Regarding the earlier lattice results, the values corre¬ 
sponding to the “pole/dipole” forms without the system¬ 
atic uncertainty estimates are displayed. Also shown are 
the QCD sum rule estimates. 


g^ = 0, = 0 and ri(0) = 72(0). Previous 

lattice results were hampered by the lack of in¬ 
formation about the heavy mass scaling behav¬ 
ior in the range of low g^’s. That generated a 
proliferation of the g^-forms to which the data 
have been fitted at large g^’s, and then extrap¬ 
olated to g^ = 0. The needed extra information 
comes from the symmetry constraints discussed in 


refs. 


allowing one to write down a parame- 


terisation similar to (^ The “pole/dipole” form 



discussed in refs, 
with symmetry constraints ]38 | 
on the new lattice results by the SPQcdR collab¬ 
oration can be found in ref. Q|. The new (still 


61 is also consistent 
More details 


® The forms used in the reanalysis of the UKQCD data in 
ref. are alail consistent with the symmetry constraints 
discussed in ISSlk^. 
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preliminary!) result for Ti(0) = T2(0) = T{0) is 

= 0.24(5) (t^) . (8) 

In fig. 1^ we compare the available results for 
the form factor T(0) computed on the lattice, 
with the values obtained by using the QCD sum 
rules . We see that the lattice estimates are 

lower than the QCD sum rule ones. Although 
the QCD sum rule method is more suitable for 
the kinematic corresponding to = 0, it is not 
clear where the discrepancy with the lattice QCD 
results comes from. More research in that direc¬ 
tion is needed (as emphasized at this conference 
in ref. (^). In addition, it should be reiterated 
that so far all the lattice studies are made in the 
quenched approximation. 

7. Conclusion 

In this talk I discussed the status of the lattice 
QCD computation of the quantities that are rel¬ 
evant to the standard determination of the shape 
of the CKM unitarity triangle (see ref. for the 
most recent update of the UTA). Quite a bit of 
progress has been made in gaining more control 
over the systematic uncertainties in the quenched 
lattice studies. A lot of work is still needed, espe¬ 
cially in unquenching the current estimates. Be¬ 
sides the “golden” quantities for the CKM trian¬ 
gle analysis, other quantities, such as the quark 
masses, decay constants, form factors, light cone 
wave functions, ... are computed on the lattice. 
Accurate determinations of these quantities is 
very useful for theoretical studies of the non- 
leptonic heavy-to-light decays. Finally, I would 
like to mention the intensive activity, within the 
lattice community, in taming the K tttt decay 
amplitudes which I did not have space nor time 
to cover in this talk. For a status report please 
see 
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